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; SUMMARY 

• A finite element was developed to adequately represent the state of stress in 

I the region around a circular hole in orthotropic material experiencing reason- 

t ably general loading. This has been achieved through a complementary virtual 

^ work formulation of the stiffness and stress matrices for a square element 

{ 

with a center circular hole. The element has been incorporated into COSMIC- 
NASTRAN as a "dunmy" element. Sample problems have been solved and these 
results are presented. 


INTRODUCTION 

Mechanical fastening is the most used joining method in present flight vehicle 
structures. This results in connection details that contain many circular holes. 
These connections are subjected to very general loading with the additional 
compllc.ttlons of load transfer through fasteners and the use of nonisotroplc 
materials in the connection. Because of the high stresses in these areas, it 
is important that an accurate evaluation of the local stresses around each 
hole be obtained to evaluate the performance of the joint. 

The finite element method is an effective way to treat problems with complex 
geometry, material, and loading. However, because of the complexity of the 
stress state around fastener holes, a large number of conventional elements 
would be required to provide an adequate representation of the hole neighbor- 
hood. This means the cost of the analysis, both in computer costs and man- 
hours, will be high. 
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A solution to this problem would be an efficient special purpose finite 
element which provides an accurate representation of the state of stress 
around a hole without the need for many conventional elements. This has 
been done previously for Isotropic materials. In 1981, under a contract 
with NASA-Langley, Lockheed-Georgla Co. developed an orthotropic hole element 
(Ref. 1). 


ELEMENT GEOMETRY 


The orthotropic hole element Is considered to be a thin plate In a state of 
plane stress with displacements on the external boundary specified so as to 
Insure displacement continuity with adjacent elements. The element Is a 
variable-noded total element and Is defined via ADUM7, CDUM7, and PDUM7 Input 
cards into NASTRAN. The user may specify the number of nodes on each of the 
three boundaries as shown in Figure 1. The nodes on the exterior of the 
clement have two cartesian degrees of freedom whereas the nodes on the hole 
wall have only a single radial degree oI freedom. The ratio of element width 
to hole diameter is determined from grid point locations and Is thus a user 
choice; however, a ratio of four Is recommended. The material properties may 
be Isotropic or orthotropic and are defined by MATl or MAT2 cards. Character- 
istics which are not defined for this element include mass properties, thermal 
loading, material nonlinearity, and differential stiffness. 

ELEMENT INPUT 

A great deal of user flexibility has been retained In the Input for this element. 
Defaults have been provided and therefore the input can be quite simple. The 
following Input is required for the element: 

o Total number of nodes 
o Element and property Identification 
o Grid point connection 
o Material Identification 

^ nnv j ^ 1. 
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The following input Is optional: 


o Number of segments on each of the three boundaries 

o Number of Integration Intervals per segment for 
each of the boundaries 

o Number of Gaussian Integration points per Interval 
for boundary integration 

o Number of Integration Intervals per volume octant 

o Number of stress terms with powers of r>- 2. 

Input for each of the four symmetry conditions. 

o Number of stress terms with powers of r< - 2. 

Input for each symmetry condition. 

o Definition of stress output locations. 

ELEMENT ADVANTAGES 

Use of the orthotropic hole element can result In a substantial savings In 
computer costs and manhours. Savings in computer costs can be obtained partic- 
ularly with models containing multiple holes because the congruency feature of 
NASTRAN can be used conveniently. This means that several hole elements of 
similar geometry can be included in the model and the stiffness matrix is only 
generated once. For complex geometries involving conventional elements, the 
congruency feature is Impractical. 

Use of the hole element can also result In substantial savings In manhours. A 
single hole element can replace dozens or even hundreds of conventional elements, 
therefore reducing the time required for modeling. Typically, an orthotropic 
hole element can be defined with eight Input cards, thus no data generators or 
preprocessors are required for these areas of the model. Another advantage 
that results In manhour savings is the stress output generated by the hole 
element. The user can obtain stress components at any location within the 
element. With a conventional model, obtaining stresses at desired locations 
can be a laborious task, because of the need to manipulate data and because of 
the large number of elements involved. 
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EXAMPLE PROBLEMS 


The exanple problemn which follow were solved with the COSMIC-NASTRAN program 
at Luckheed-Georgla Company with the orthotropic hole element implemented. 
They are Intended to demonstrate the performance of the element In a variety 
of applications. 


Isotropic Plate with Uniaxial Tension 

The orthotropic hole element can be used in isotropic applications by simply 
specifying isotropic elastic coefficients. This is demonstrated by presenting 
the solution for an infinite isotropic plate with unit uniaxial tension at 
infinity (see Figure 2). The finite element model is shown in Figure 3. The 
plate width is 20 times the hole diameter, which represents an essentially 
infinite plate. Poisson's ratio is set at 0.3205. 

Table 1 presents results for the three stress components at four values of 
radius and four values of 6. The hole element stiffness and stress matrices 
were based on a geometry of 8 exterior nodes and 12 interior nodes (on the 
hole boundary). The hole boundary is given by r = 1. Note the excellent sat- 
isfaction of the zero radial stress boundary condition; the results were typi- 
cally 10 times the magnitude of external loading. The circumferential stress 
results are also very good, with an error of only 0.2% on the maximum tension 
stress concentration factor and 1.5% on the maximum compression stress. The 
zero shear stress condition is perfectly satisfied since it was included in 
the original stress assumptions. 

The stresses away from the hole boundary are also very well represented by 
the hole element results. 
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Isotropic Plate with Half Cosine Loaded Hole 


A load transfer problem for which a comparison solution Is available (Ref . Is 
that of an Infinite Isotropic plate with Poisson's ratio equal to 0.25, loaded 
with uniform traction on one end and a cosine variation of radial tractions on 
the hole boundary (see Figure 4). The total load on the end Is equal to 2P, 
and the radial tractions on the loaded hole boundary are given by 

T “ 4P cos <b t w 

t — 2< '*' ■= 2 

The fln-ice element model Is shown in Figure 3. Table II presents results for 
«-Se ,.eglcn around the hole, ^a-ied on a geometry of 8 exterior nodes and 16 
Interior nodes. The value of P Is set equal to 1.0. 

Stresses at the hole boundary are very well represented by the finite element 

results. The radial stresses follow the cosine curve over the loaded region 

-2 

and are of the order of 10 on the unloaded portion of the boundary. The 
variation of circumferential stress Is reasonably accurate, with an error of 
3Z on the maximum str.:ss concentration. 

In the region around the hole, the results from the hole element solution 
continue to give a good description of the stress state. 

^5° Laminated Plate 

The example problems with Isotropic material serve to confirm the basic theoret- 
ical approach and to partially validate the computer solution. There remains to 
demonstrate satisfactory performance In applications with orthotropic materials. 
The first material selected has material properties as follows: 

E - E - E, G - 1.697528E, v = v =■ 0.735 
X y ’ xy yx xy 

Note that even though E “ E , this material is rather "far" from isotropic, 

X y 

with some unspecified measure, due to the large value of Poisson's ratio and 
the fact that E, G, and v do not satisfy the Isotropic relationship. These 
properties describe the ^5° graphlte/epoxy laminate considered In Reference 4. 
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The first loading is unit tensile tractions in the y-direction, applied as 
shown in Figure 2 with the finite element model given in Figure 3. The 
computed radial stresses at the hole boundary were of the order 10 ^ times 
the exterior loading thus satisfying traction boundary conditions. Figure 5 
shows the variation of circumferential stress around the hole boundary for a 
hole element geometry consisting of 8 exterior nodes and 12 interior nodes. 

Also shown is the exact solution for an infinite plate, given by 

''’e “ 3.05786 cos^ 9-1 

2.8809072(cos^ 6 - cos^ 0) + 1 

which can be deduced from Reference 5, p. 175. The orthotropic hole element 
results contain Che essential feature that the maximum stress concentration 
does not occur at the familiar location of 0 = 0°. The error in maximum stress 
is only 0.7%. 

The hole loaded with half-cosine distributed tractions has also been solved; 
see Figure 4. The resulting radial stresses on the boundary are essentially 
Identical to Che results achieved in the case of isotropic material. The 
maximum circumferential stress on the hole boundary is Tg * 0.79 at 125°, 
which can be ompared with the hole element results for isotropic material of 
Tq = 0.84 at 0= 85°. Again, it is seen that the location of maximum stress 
has shifted around the circle from the familiar location. 

Plywood 

The next material selected has properties as follows: 

E = E. E - 2E, G - 0.11667E, v - 0.036 
y X ’ xy xy 

These properties describe a pl}wood material for which Reference 5 contains many 
results for an infinite plate. Note that the plywood properties ace quite 
different from the ;H5° lamination properties, both in magnitudes and relation- 
ships, This provides further checks on the ability of Che hole element to 
handle different types of orthotroplc material. 
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The types of loading considered are unit tensile tractions in the y-direction 
(Figure 2), unit compressive tractions in the x-direction with constraint in 
the y-direction (Figure 6a) , and uniform unit internal pressure on the hole 
boundary (Figure 6b) . For all cases, the finite element model is shown in 
Figure 3; and the element geometry contains 8 exterior nodes and 24 interior 
nodes . 

For the uniaxial tension, there is the usual excellent satisfaction of traction 
boundary condition, with of the order 10 Circumferential stresses around 
the hole boundary are shown in Figure 7. Also shown are the exact results for 
an infinite plate, given by 

Tq = 4.8557 cos^ e - 0.7071 

-6.9994 cos^ 0 + 7.4994 cos^ 6+0.5 

which can be deduced from Reference 5 p. 175. The general features of the variation 
of stress are demonstrated by the hole element results. 

The constrained compression results satisfy the hole boundary conditions, with 
Tp of the order 10 Circumferential stresses around the boundary are shown 
in Figure 8. The exact results tor infinite plate are given by 

r. = 1.1155 - 6.5174 sin^e 

V 4 2 

-13.9989 sin 0 + 12.9989 sin 0+2 

deduced from Reference 5 p. 179. 

The internal pressure results provided boundary radial stresses of the magnitude 
-1.00002, which is again excellent satisfaction of boundary conditions. Circum- 
ferential stresses around the boundary are shown in Figure 9, compared to 
infinite plate exact results given by Reference 5 p. 173. 

= 4.8829 - 15.8432 sin^ 6 + 13.9989 sin^6 

^ 2 4 

2 + 12.9989 sin 0 - 13.9989 sin 0 




Orthotropic Laminates with Pin-Loaded Holes 


A class of problems which is of particular interest to the aircraft industry 
is that of pin-loaded holes in orthotropic laminates. Finite element solutions 
for a variety of graphite epoxy laminates were presented by Crews (Ref. 6). The 
finite element model used for these analyses was quite detailed, as shown in 
Figure 10. These results are used here as a basis of comparison to results 
obtained using the orthotropic hole element. 

A NASTRAN model was constructed as shown in Figure 11. The particular case 
presented here is for W/d * 20. The load is applied at the center of the 
fastener and is reacted on both edges of the model in order to more closely 
approximate an infinite sheet. A feature of Lockheed-Georgia' s version of 
COSMIC-NASTRAN was Invoked to properly account for the contact area of the 
fastener. The hole element geometry consisted of 16 exterior nodes and 20 
interior nodes. Six different laminates were considered. These were (1) quasi- 
isotropic (2) 0*^ (3) 90° (4) [o°/90 ]s (5) [+45°]s and (6) [o°/+45°]s. 

Figures 12-17 present plots of the tangential and radial stress components 
obtained by the hole element for these respective materials in comparison to 
the results obtained by Crews using a detailed classical model. 

These plots indicate that a reasonable correlation was obtained for the tangential 
stress distribution at the hole boundary. Peak tangential stress differed by no 
more than 12%. Excellent correlation of radial stress distribution was obtained 
for all cases. 


Multiple Holes 

Figure 18 shows a finite element model of a plate containing five holes. For 
economy in modeling, it is desirable to develop the hole element stiffness 
matrix only once and then transform this matrix as necessary to represent all 
the hole elements. This is possible in NASTRAN by using the congruency feature, 
and this feature has been checked by determining the solution to Figure 18 with 
tension applied parallel to the row of holes. The material is Isotropic with 
Poisson's ratio of 0.3205. 
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The maximum stress concentration factor Is given In Reference 7, p. 194 for an 
Intermediate hole In an Infinite row of holes. In order to simulate the 
Infinite row, the displacements of each loaded end of the model were constrained 
to be uniform across the width of the plate. The finite element result for 
stress concentration Is 2.78, based on a hole element geometry containing 8 
exterior nodes and 12 Interior nodes. Reference 7 gives 2.65; there Is a 5% 
difference In results. 


( 

I 
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Concluding Remarks 

The performance of the orthotropic hole element in solving a variety of problems 

has demonstrated that it can produce acceptable results with significant advantages 

In computer costs, model preparation time, and post-processing of data. Even i 

though the emphasis In this paper has been on problems concerning fastener holes, 

the use of this element Is applicable to any Isotropic or orthotropic plate , 

containing circular holes. 

I 
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Figure 2 


Unit Tensile Traction in the Y-Dlrection 







hole element 
results 



Figure 5. Circumferential Stress at the Open Hole Boundary for 
a +45° L-'minate with Unit Uniaxial Tensile Traction 







exact solution for 
infinite plate 



Figure 7. Circumferential Stress at the Open Hole Boundary 
for Plywood with Unit Uniaxial Tensile Traction 




Figure 8. Circumferential Stress at the Open Hole Boundary 
for Plywood with Unit Constrained Compression 
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Figure 10. Detailed Finite Element Model for Analysis 

of Orthotropic Laminates with Pin-Loaded Holes 
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Figure 12. Stresses in the Quasi-Isotroplc Laminate 
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Figure 18. Finite Element Model for a Plate 
with Multiple Holes 
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